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ABSTRACT 

Bolometric luminosities and Eddington ratios of both X-ray selected broad-line (Type-1) and 
narrow-line (Type-2) AGN from the XMM-Newton survey in the COSMOS field are pre- 
sented. The sample is composed by 929 AGN (382 Type-1 AGN and 547 Type-2 AGN) and it 
covers a wide range of redshifts, X-ray luminosities and absorbing column densities. About 
65% of the sources are spectroscopically identified as either Type-1 or Type-2 AGN (83% 
and 52% respectively), while accurate photometric redshifts are available for the rest of the 
sample. The study of such a large sample of X-ray selected AGN with a high quality multi- 
wavelength coverage from the far-infrared (now with the inclusion of Herschel data at 100 /urn 
and 160 fim) to the optical-UV allows us to obtain accurate estimates of bolometric luminosi- 
ties, bolometric corrections and Eddington ratios. The fct,oi-£boi relations derived in the present 
work are calibrated for the first time against a sizable AGN sample, and rely on observed red- 
shifts, X-ray luminosities and column density distributions. We find that k\, \ is significantly 
lower at high Lboi with respect to previous estimates by Marconi et al. (2004) and Hopkins et 
al. (2007). Black hole masses and Eddington ratios are available for 170 Type-1 AGN, while 
black hole masses for Type-2 AGN are computed for 48 1 objects using the black hole mass- 
stellar mass relation and the morphological information. We confirm a trend between k\, \ and 
^Edd, with lower hard X-ray bolometric corrections at lower Eddi ngton ratios for both Type- 1 
and Type-2 AGN. We find that, on average, Eddington ratio increases with redshift for all 
Types of AGN at any given Mbh, while no clear evolution with redshift is seen at any given 

£bol- 
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1 INTRODUCTION 

It is widely accepted that the central engine of AGN are accreting 
supermassive black holes (SMB Hs) at the center of galaxies wit h 
masses of the order of 10 6 ~ 9 M o dSalpetedl964l : lLvnden-Bellll969h . 
Locally, the SMB H mass correlates with the mass of the bulg e 
of the host-galaxy dMagorrian et al.ll 19981 ; IMarconi & HuntH2003l). 
with the velocity dispersi on of the bulge ( Ferrarese & Merrittl 
120001: iTremaine et al.|[2 002). and with the luminosity of the bulge 
dKormendv & Richstondll995r) . The existence of these correla- 
tions implies that the growth of the SMBH is tightly linked with 
the galaxy evolution, playing a crucial role in the star-formation 
history of the galaxy itself. The feedback between SMBH and 
host-galaxy is therefore a pivotal ingredient that has to be taken 
into account in SMBH/galaxy format i on and co-evolution studies 
(see ISilk & Reesl 1 19981 : iFabianl 1 19991 : Iciotti & van Albadal l200ll : 
Di Ma tteo et alJl2005D ^ A fundamental question is how the AGN 
energy detected as radiation is produced. AGN are a broad-band 
phenomenon, hence a multi-wavelength approach is mandatory in 
order to understand the physics that underlie the AGN emission. 
Several studies have been performed on the shape of the SED, as 
parameterized by the correlation between the a ox index, defined 
as q ox = -Log [L 2 kev/ L 7 srm a]/2.605, a nd th e optical luminosity 
dTananbaum et al J| 197 9 : Zamora ni et al.ll 1 98 ll; I Vignali et al.ll2003l: 
Steffe n et al.ll2006l: I Just et al.ll2007l : I Young et alJl2009l : lLusso et al.l 
20ld:lMarchese et alJl2012[T or between a m and t he Eddington ra 



tio ( Vasudevan & Fabian 12007|; iKellv et alj|2008l: Vasudevan et al. 
120091 : IVasudevan & Fabianl 120091 : iTrump et all 1201 ll) . How o~ 
evolves with luminosity may provide a first hint about the nature of 
the dominant energy generation mechanism in AGN. It is also a first 
step towards an estimate of the AGN bolo metric luminosity func- 
tion dHopkins et al]|200l here after H07:lshankar et al.ll2009h and 
the m ass function of SMBHs dShankar et alj|2004l : IMarconi et al.l 
1 20041 hereafter M04). All these works consider a model intrinsic 
quasar SED, which is described with a series of broken power- 
laws, similar to those objects in bright optically selected samples 
dElvis et al.|[l994l : iRichards et al]|2006h . M04 and H07 have also 
studied the relationship between the bolometric correction, fcboi, as a 
function of the bolometric luminosity, L bo i, in different bands (e.g., 
the B-band at 0.44pm, the soft and the hard X-ray bands at [0.5-2] 
keV and [2-10] keV, respectively). Therefore, it is of fundamental 
importance to verify these correlations considering statistically rel- 
evant samples of both broad-line (Type-1) and narrow-line (Type- 
AGN over a wide range of redshift and luminosities. 
We analyze the dependence of k bo] on L bo] in the B-band, and 
in the soft and hard X-ray bands using a large X-ray selected 
sample of both Type-1 and Type-2 AGN in the Cosmic Evolu- 
tion Survey (COS MOS) field dScoyille et al.l2007h from the XMM- 
COSMOS sample dHasinger et aT]|2007h . The COSMOS field is a 
unique area for its deep and wide multi- wavelength coverage: radio 
with the VLA, infrared with Spitzer and Herschel, optical bands 
with Hubble, Subaru, SDSS and other ground-based telescopes, 
near- and far-ultraviolet bands with the Galaxy Evolution Explorer 
(GALEX) and X-rays with XMM-Newton and Chandra. The spec- 
troscopic coverage with VIMOS/VLT and IMACS/Magellan, cou- 
pled w ith the reliable photo metric redshifts derived from multiband 
fitting dSalvato et al.ll201ll herefter SI 1), allows us to build a large 
and homogeneous sample of AGN with a well samp led spectral 
coverage and to keep selection effects under control. In lLusso et al.l 



d2O10L hereafter L10) bolometric corrections and Eddington ratios 
for hard X-ray selected Type- 1 AGN in the COSMOS field are pre- 
sented. That work showed that the bolometric parameters are useful 
to give an indication of the accretion rate onto SMBH (i.e., a lower 
bolometric correction corresponds to a lower Ed d ington ratio, and 
vice- v ersa; see also I Vasudevan & Fabianl [2009: Vas udevan et all 
2009, hereafter V09a and V09b respectively). However, in L10, 
V09a and V09b the study was mainly focused on the bolometric 
output of Type-1 AGN, while an alternative approach is needed 
in order to study the bolometric parame ters and Eddington ratios 
of Type-2 AGN. I Vasudevan et all feoid) (hereafter V10) explore a 
method for characterizing the bolometric output of both obscured 
and unobscured AGN by adding nuclear IR and ha rd X-ray lumi- 
nosities (as originally proposed bv lPozzi et al .l2007t) . They also es- 
timate Eddington ratios using black hole mass estimates from the 
black hole mass-host galaxy bulge luminosity relation for obscured 
and unobscured AGN finding a significant minorit y of higher accre- 
tion rate objects amongst high-absorption AGN. Lusso et al. I d201lh 
(hereafter Lll) present an SED-fitting method for characterizing 
the bolometric output of obscured AGN, and giving also an esti- 
mate of stellar masses and star formation rates for Type-2 AGN. 
In this analysis the mid-infrared emission is utilized as a proxy 
to constrain the thermal emission, and it is used, in conjunction 
with the X-ray emission, to give an estimate of the bolometric lu- 
minosity. Taking advantage of these results, we are able to build 
up a homogeneous analysis to further study the bolometric out- 
put of X-ray selected AGN at different absorption levels. More- 
over, several improvements are included in this study with respect 
to L10 and L ll. Firstly, the in clusion of Herschel data at 100 pm 
and 160 pm dLutz et alj|201lh in the SED-fitting code is of fun- 
damental importance to better constrain the far-infrared emission 
and, therefore, the AGN emission in the mid-infrared for Type-2 
AGN. Photometric redshifts are taken into account in order to ex- 
tend our sample also to fainter magnitudes using the updated re- 
lease of photometric redshifts provided by SI 1. Finally, we have 
included the H band photometry now available in the COSMOS 
field; this allows us to further increase the coverage in the near- 
infrarecj|. Black hole masses and Eddin gton ratios for Type-2 AGN 
are estimated through scaling relations (Haring & Rix 2004) using 
a Monte Carlo method in order to account for uncertainties in stellar 
masses, bolometric luminosities, as well as the intrinsic scatter in 
the M B h - M, relation. The morphological information is also taken 
into acco unt in order to estima te the bulge-to-total luminosity/flux 
ratio (see lSimmons et al]|201 ll and Sect.|5]of the present paper for 
further details) which is used to further constrain the black hole 
mass estimate. 

This paper is organized as follows. In Sect. [2] we report the 
selection criteria for the hard X-ray selected samples of Type-1 
and Type-2 AGN used in this work. Section [3] presents the multi- 
wavelength data-set, while Section [4] describes the methods used 
to compute intrinsic bolometric luminosites and bolometric correc- 
tions for Type-1 and Type-2 AGN. Section |5]describes black hole 
mass estimates and Eddington ratios, while in Sect. [6] we discuss 
our findings. In Sect.|7]we summarize the most important results. 

We adopted a fiat model of the universe with a Hubble 
constant Hp = 70 k m s -1 Mpc~', Q. M = 0.27, Q A = 1 - fl M 
dKomatsu et alj2009l) . 



1 In the following we will adopt the nomenclature "Type-2" AGN, "not 2 These data were re duced using similar methods as described in 

Type-1", or "not broad-line" AGN for the same population of sources. iMcCracken et alJfcOlOl) . 
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Table 1. Summary of the photometric classification criteria. 



Sample 


Total 


SED-Type a 


N b 


Photo-z 


330 


Galaxy-dominated SED 
AGN-dominated SED 


263 
67 


Type-1 Spectro-z 


315 


(ill r\ v \;_,rirMTi i n n tf^ri v! H Ti 

AGN-dominated SED 
no photo-z Class 


3Q (\0WA 

271 

5 


Type-2 Spectro-z 


284 


Galaxy-dominated SED 
AGN-dominated SED 
no photo-z Class 


257 

25 (9%) 
2 



[a] Photometric classification according to S 1 1 . Templates 19-30 have 
AGN-dominated SEDs, while templates 1-18 and 101-130 have 
galaxy-dominated SEDs. 

[b] Number of sources in each sub-sample: galaxy-dominated SEDs, 
AGN-dominated SEDs, and sources without photometric classification. 

2 THE DATA SET 

The Type-1 and Type-2 AGN samples discussed in this paper are 
extracted from the XMM-COSMOS catalog which comprises 1822 
point-like X -ray sources detected by XMM-Newton ove r an area 
of ~ 2 deg 2 faasinger et al]|2007l: Icappelluti et alj|2009b. A ll the 
details about the catalog are reported in iBrusa et al.l (120101) . We 
consider in this analysis 1577 X-ray selected sources for which 
a reliable optical counterpart can be associated (see discussion in 
IBrusa et al.l |2010L Table lfl A not negligible fraction of AGN 
may be present among optically faint X-ray sources without op- 
tical spectroscopy; therefore, in order to extend both Type-1 and 
Type-2 AGN samples to fainter luminosities, we have considered 
the updated release of the photometric catalog provided by Sll. 
Photometric classification from 1 to 30 include galaxy-dominated 
SEDs (early type, late type and ULIRG galaxies), low-and-high lu- 
minosity AGN SEDs and hybrids created assuming a varying ratio 
betwe en AGN and galaxy templates (see Table 2 in lSalvato et al.l 
2009 for details). Sources coded from 100 to 130 are reproduced 
by host-galaxy domina ted SED (elliptica l, spiral and star-forming 
galaxies, see Fig. 1 in lllbert et al.ll2009h . We differentiate photo- 
metric classification using the best-fit template, separating those 
dominated by AGN emission and galaxy emission. Selection cri- 
teria for the Type-1 and the Type-2 AGN samples are described in 
the following Section. 

2.1 Type-1 and Type-2 AGN samples 

We have selected 971 X-ray sources detected in the [2-10] keV 
band at a flux larger than 3 x 10~ 15 erg s~'crrT 2 (see IBrusa et al.l 
l2O10t) . From this sample, 315 objects are spectroscopically clas- 
sified broad-line AGN on the basis of broad emi ssion lines 
(FWHM > 2000 km sec -1 ) in their optical spectra (see lLillv et al.l 
l2007t iTrump et all l2009bh . We will refer to this sample as the 

3 The multi-wavelength XMM-COSMOS catalog can be retrieved 
from: http://www.mpe.mpg.de/XMMCosmos/xmm53_release/ version 1 st 
November 2011. 



"spectro-z" Type-1 AGN sample. We consider Type-2 AGN all 
the rest, including Seyfert 2 AGN, emission-line galaxies, and 
absorption-line galaxies. The remainder 284 AGN with optical 
spectroscopy do not show broad emission lines (FWHM< 2000 
km s~') in their optical spectra: 254 are objects with either unre- 
solved, high-ionization emission lines, exhibiting line ratios indi- 
cating AGN activity, or not detected high-ionization lines, where 
the observed spectral range does not allow to construct line diag- 
nostics. Thirty are classified absorption-line galaxies, i.e., sources 
consistent with a typical galaxy spectrum showing only absorption 
lines. We will refer to this sample as the "spectro-z" Type-2 AGN 
sample. 

In order to extend our Type-1 and Type-2 AGN samples to 
fainter magnitudes, we proceed as follow. We have selected all 
sources with a best-fit photometric classification consistent with an 
AGN-dominated SED (i.e., 19 < SED - Type < 30; see Sll, for 
details). In the following, we assume that 67 X-ray sources, clas- 
sified by the SED fitting with an AGN-dominated SED are Type- 1 
AGN. We will refer to this sample as the "photo-z" Type-1 AGN 
sample. From the total photo-z sample we additionally selected all 
AGN with a best-fit photometric classification inconsistent with a 
broad-line AGN SED. We will refer to this sample (i.e., 263 X-ray 
sources with SED - Type > 100 & SED - Type < 19, see TableQ} 
as the "photo-z" Type-2 AGN sample. 

The final Type-1 AGN sample used in our analysis comprises 
382 X-ray selected AGN (315 from the spectro-z sample and 67 
from the photo-z sample), while the final Type-2 AGN sample com- 
prises 547 X-ray selected AGN (284 from the spectro-z sample and 
263 from the photo-z sample). These samples span a wide range of 
redshifts and X-ray luminosities. Only 12 Type-2 AGN (2%, 9 and 
3 AGN from the spectro-z and photo-z sample, respectively) have 
Log Lp-iojkev < 42 erg s . Sources with Lp-iojkev < 42 might be 
interpreted as star-forming galaxies, however the number of these 
objects is very small and their inclusion in the Type-2 class does 
not affect any of the results. 

We have estimated the expected contamination and incom- 
pleteness in the classification method, on the total Type-1 and Type- 
2 AGN samples, by checking the distribution of the photometric 
classification for the spectroscopically identified Type-1 and Type- 
2 AGN samples. The large majority of the broad emission line AGN 
in the spectro-z sample are classified as Type-1 AGN by the SED 
fitting as well (271/315; 86%), while the number of spectroscopic 
Type-1 AGN which have SED-Type > 100 or SED-Type < 19 is 
relatively small (39 sources). Similarly, a good agreement between 
the two classifications is present also for the Type-2 AGN; for 91% 
(257/284) of the 284 Type-2 sources with spectroscopic redshift, 
the spectroscopic classification is in agreement with the SED-based 
classification. If we make the reasonable assumption that the frac- 
tion of agreement between the two classifications is the same also 
for the sample for which we do not have a spectroscopic classifica- 
tion, we find that our total Type-1 AGN sample is expected to be 
contaminated (i.e. Type-2 AGN misclassified as Type 1 AGN from 
the SED analysis) at the level of ~1.6% (6 objects) and incomplete 
(i.e. Type-1 AGN misclassified as Type-2 AGN, and therefore not 
included in our Type-1 sample) at the level of ~9.2%. The same 
fractions for the total Type-2 sample are ~6.4% (contamination) 
and ~1.1% (incompleteness). In Table [T] a summary of the selec- 
tion criteria and the numbers of AGN in the various classes are 
presented. 

The redshift distributions of the total, spectroscopic and pho- 
tometric Type- 1 and Type-2 AGN samples are presented in Figure[T] 
(left panel). The median redshift of the total Type-1 AGN sample 
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1 2 3 1 2 3 4 

redshift redshift 



Figure 1. Redshift distribution of the hard X-ray selected Type-1 (left panel) and Type-2 (right panel) AGN samples considered in this work. The hatched 
histogram shows the redshift distribution for the sample of spectroscopically identified sources, while the filled histogram is the redshift distribution for the 
sources without spectroscopic redshift. The total sample is reported with the open histogram. 



is 1.51 (the mean redshift is 1.57, with a dispersion of 0.70). The 
median redshift of the spectro-z sample is 1.45, while the median 
redshift of the photo-z sample is 1.92. The median redshift of the 
total Type-2 AGN sample is 0.96 (the mean redshift is 1.10, with a 
dispersion of 0.64). The median redshift of the spectro-z sample is 
0.81, while the median redshift of the photo-z sample is 1.41. 



3 SPECTRAL ENERGY DISTRIBUTIONS 

We have collected the multi-wavelength information from mid- 
infrared to hard X-rays as in L10 and Lll. The observations in 
the various bands are not simultaneous, as they span a time in- 
terval of about 5 years: 2001 (SDSS), 2004 (Subaru and CFHT) 
and 2006 (IRAC). Variability for absorbed sources is likely to be a 
negligible effect, which is probably not the case for Type- 1 AGN. 
Therefore, in order to reduce variability effects, we have selected 
the bands closest in time to the IRAC observations (i.e., we ex- 
cluded SDSS data, that in any case are less deep than other data 
available in similar bands). All the data for the SED computa- 
tion were shifted to the rest frame, so that no K-corrections were 
needed. Galactic reddening has been taken into account: we used 
the selective attenuation of the stellar continuum k(X) taken from 
Table 1 1 of | Capak et all d2007t) . Galactic extinction is estimated 
from lSchlegel et al.ldl998l) for each object. We decided to consider 
the near-UV GALEX band for Type-1 and Type-2 AGN with red- 
shift lower than 1, and far-UV GALEX band for sources with red- 
shift lower than 0.3 in order to avoid Lya absorption from fore- 
ground structures. In th e far-infrared the inclusion of Herschel data 
at 100 /mi and 160 yum dLutz et alj|201ll) better constrain the AGN 
emission in the mid-infrared. The number of detections at 100 yum 
is 63 (16%, 59 spectro-z and 4 photo-z) for the Type-1 AGN sam- 
ple, while is 98 (18%, 73 spectra-z and 25 photo-z) for the Type-2 
AGN sample. At 160 [im the number of detections for the Type- 
1 AGN sample is 56 (15%, 52 spectro-z and 4 photo-z), while is 
87 (16%, 63 spectra-z and 24 photo-z) for the Type-2 AGN sam- 
ple. Count rates in the 0.5-2 keV and 2-10 keV are converted into 
monochromatic X-ray fluxes in the observed frame at 1 and 4 keV, 
resp ectively, using a Galactic column density Nh = 2.5 x 10 20 cm -2 
(see lDickev & LockmaiJl990l : lKalberla et alj2005T) . We have com- 
puted the integrated unabsorbed luminosity in the [0.5-2]keV and 
[2-10]keV bands for both Type-1 and Type-2 AGN samples. For a 
sub-sample of 100 Type-1 AGN (26%) and 240 Type-2 AGN (44%) 



we have an estimate of the co l umn d ensity Nh from spectral anal- 
ysis (see iMainieri et alj|2007l , l20ld) . while for 282 Type-1 AGN 
and 307 Type-2 AGN abso rption is estimated from hardness ratios 
(HR; see lBrusa et alj20ich . The integrated intrinsic unabsorbed lu- 
minosity is computed assuming a power-law spectrum with slope 
r = 2 and T = 1.7 for the [0.5 -2]keV and [2-10]keV bands, respec- 
tively dCappelluti et alj2009h . 

In Figure [2] we show the distribution of column densities for 
the Type-1 and Type-2 AGN samples. Two-hundred and twenty- 
four Type-1 AGN (~ 64%) and 87 Type-2 AGN (~ 16%) do not 
require absorption in addition to the Galactic one. The mean N H 
value is 7.4x 10 20 citT 2 for the Type-1 AGN and ~ 10 22 cirT 2 for the 
Type-2 AGN. Sixty percent (322/547) of the Type-2 AGN sample 
and 10% (39/382) of the Type-1 AGN sample have Log N H > 22 
cnT 2 . The average shift induced by the correction for absorp- 
tion in the Type-1 sample is, as expected, small in the soft band 
(ALog L[o.5-2]kev) = 0.10 ± 0.01, negligible in the hard band.The 
same shift in the Type-2 sample is (ALog L[o.s-2jkev) = 0.30 ± 0.04 
in the soft band, while it is (ALog L [2 -io]keV> = 0.033 ± 0.004 for 
the hard band. 

We have computed the individual rest-frame SEDs for all 
sources in the sample, following the same approach as in L10. 
For the computation of the bolometric luminosity for Type- 1 AGN 
we need to extrapolate the UV data to X-ray gap and at high X- 
ray energies. We have extrapolated the SED up to 1200 A with the 
slope computed using the last two rest-frame optical luminosity 
data points at the highest frequency in each SED (only when the 
last optical-UV rest frame data point is at A > 1200A). Then, we 
assume a power law spectrum to 500A, as measured by HST obser- 
vations for radio-quiet AGN (f v cc y~' 8 , see lZheng et all 19971) . We 
then linearly connect (in the log space) the UV luminosity at 500 A 
to the luminosity corresponding to the frequency of 1 keV. 



4 BOLOMETRIC LUMINOSITIES 

The mid-infrared luminosity is considered an indirect probe of 
the accretion disk optical/U V luminosity (see IPozzi et all 1 20071 , 
l201fi IVasudevan et alj[2O10h . The nuclear bolometric luminosity 
for Type-2 AGN is then estimated by using the same approach as 
in Lll, whereas the sum of the infrared and X-ray luminosities 
are used as a proxy for the intrinsic nuclear luminosity (L bo] = 
Lir + Lx). The main purpose of the SED-fitting code is to disentan- 
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Log N H [cm -3 ] Log N H [cm -2 ] 

Figure 2. Column density distribution of the Type-1 (left panel) and Type-2 (right panel) AGN samples (black open histogram). About 67% of the Type-1 
AGN sample and 16% of the Type-2 AGN sample have N\i values consistent with the Galactic one. Note that the first bin of in the left panel is much higher 
than what is plotted in the histogram (206 spectro-z and 40 photo-z have Nh = 20.5 cm -2 ), while for the right panel 56 spectro-z and 33 photo-z have Nu 
consistent with the Galactic value. The hatched histogram shows the distribution for the sample of spectroscopically identified sources, while the filled 
histogram is the Nu distribution for the sources without spectroscopic redshift. Sixty percent (322/547) of the Type-2 AGN sample and 10% (39/382) of the 
Type-1 AGN sample have Log Nu ^ 22 cirT 2 . 
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Figure 3. Examples of SED decompositions. Black circles are the observed photometry in the rest-frame (from the far-infrared to the optical-UV). The blue 
long-dashed, black solid and dotted lines correspond respectively to the starburst, AGN and host-galaxy templates found as the best fit solution. The red line 
represents the best-fit SED. The top three panels are Type-1 AGN, while the bottom three are Type-2 AGN. 



gle the various contributions (star-burst, AGN, host-galaxy emis- 
sion) in the observed SEDs by using a standard x 1 minimization 
proce dure. The code is base d on a large set of star-bu rst templates 
from lCharv & Elbaj j200lT) and bale & Heloul d2002h . and galaxy 
templates from the lBruzual & Cha riot (2003) code for s pectral syn- 
thesis models, while AGN templates are taken from Isilva et all 
d2004l) . These templates represent a wide range of SED shapes and 
luminosities and are widely used in the literature. After perform- 
ing the SED-fitting, only the nuclear component of the best-fit is 



integrated. Hence, the total IR luminosity Lir is obtained integrat- 
ing the nuclear template between 1 and 1000 jim. To convert this 
IR luminosity into the nuclear accretion disk luminosity, we ap- 
plied the correction factors to account for the torus geometry and 
the anisotropy (~ 1.7, see IPozzi et afll2007l I20T0I) . The total X- 
ray luminosity Lx is estimated by integrating the X-ray SED in the 
0.5-100 keV range. 

The photometric data used in the SED-fitting code, from 
low to high frequency, are: Herschel-V 'ACS bands (160 j-im and 
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100 fan), MlPS/Spitzer (24 yum and 70 lira), 4 IRAC bands (8.0 yi/m, 
5.8 yum, 4.5 yum and 3.6 yum), K s CFHT, # CFHT, 7 UKIRT, op- 
tical broad-band Subaru (Bj, V, g + , r + and z + ) and CFHT («'*, «*) 
bands. The starburst component is used only when the source is de- 
tected at wavelengths longer than 24 fim rest-frame. Otherwise, a 
two components SED-fit is used. Eighteen is the maximum number 
of bands adopted in the SED-fitting (only detection are considered). 
Figure[3]shows the multi-wavelength photometry and the model fits 
for six AGN, three Type-1 AGN in the top panels and three Type-2 
AGN in the bottom panels. The three components adopted in the 
SED-fitting code, starburst, AGN torus and host-galaxy templates, 
are shown as a blue long-dashed line, black solid line and dotted 
line, respectively. 

The bolometric luminosity for Type-1 AGN is usually com- 
puted by integrating the observed SED, as described in Sect. [3] in 
the Log v - Log (vL y ) rest-frame plane from lyum to 200 keV. The 
choice to neglect the infrared bump is motivated by the fact that 
nearly all photons emitted at these wavelengths by the AGN are re- 
processed optical/UV/soft X-ray photons; in this way we avoid to 
count twice the emission reprocessed by dust (see M04). However, 
given that we want to compare bolometric parameters for Type- 
1 and Type-2 AGN, we have decided to use the SED-fitting code 
above described to compute bolometric luminosities for both sam- 
ples. 

In order to compute the bolometric correction we used the 
standard definition 



30 



3 

10 



373' Type-1 'AGN 
488 Type-2 AGN, 




44 46 

Log L bol [erg s" 1 ] 

Figure 4. Histogram of bolometric luminosities for Type- 1 AGN (blue filled 
histogram), and for Type-2 AGN (red open histogram) with AGN best-fit 
available from the SED-fitting code. The two dashed lines show the median 
Lbol values for Type-1 AGN L bo i = 3.2x 10 45 erg s -1 (right), and for Type-2 
AGN L bol = 7.1 X 10 44 erg s" 1 (left). 



(1) 



where L ban d ' s me luminosity in soft and hard bands and in the B- 
band at 0.44/jm. The luminosity in the B-band is computed only 
for Type-1 AGN, since for Type-2 AGN the emission in the optical 
is mainly from the host-galaxy. Bolometric luminosities for Type- 
1 and Type-2 AGN samples are reported in Figure |4] The mean 
bolometric luminosity value for Type-1 AGN is (Log L bol ) = 45.50 
with a dispersion of 0.58, while for Type-2 AGN (Log L bo i) = 44.85 
with a dispersion of 0.70. Using the SED-fitting approach we found 
that 9 Type-1 AGN are best-fitted either with only galaxy tem- 
plate (7 objects mostly at redshift higher than 2), or with galaxy 
plus star-burst template (2 objects, all of them have Herschel data 
at 100/160 /jm). Therefore, bolometric AGN luminosities are not 
available for these sources. For 59 Type-2 AGN (10% of the main 
sample) the SED-fitting code is not able to fit the mid-IR part 
of the SED with the AGN component. Twenty-two of these 59 
Type-2 AGN are well fitted only with a galaxy component. This 
sample is predominantly at high redshift, highly obscured with 
N H > 10 23 cirr 2 and none of these sources have a detection in 
the far-IR, most of them neither at 24 /jm band. Given that they 
are mainly at high redshift, all bands are shifted towards high fre- 
quencies, therefore the torus component is not needed. For the re- 
maining 37 Type-2 AGN the best-fit is composed by the galaxy 
component in the optical and the star-burst component in the far- 
IR. All these sources present MIPS detection at 24/70 yum and/or 
Herschel data at 100/160 yum. These Type-2 AGN are at relatively 
moderate redshift (z < 1.5) and X-ray luminosities in the range 
40.5 < Lp-iojkev ^ 44. Summarizing, the SED-fitting for 2% of the 
Type-1 AGN sample and for 1 1% of the Type-2 AGN sample is not 
able to recover the AGN component in the mid-infrared. 

We have estimated the uncertainties on bolometric luminosi- 
ties comparing the bolometric luminosities computed by the SED- 
fitting code and those obtained integrating the rest-frame SED from 
lyum to the X-ray. Obviously, this test can be only applied to the 



Type-1 AGN sample. In Figure[5]the comparison between L bo i from 
the integrated SED and from the SED-fitting is presented. We find 
that the bolometric luminosities are scattered along the one-to-one 
correlation with a 1-sigma scatter of ~ 0.24 dex (after performing 
a 3.5cr clipping, 4 objects have been removed). Assuming that the 
errors associated to the two different ways to compute L bol are of 
the same order of magnitude, we can estimate the 1 a uncertainty 
on L bo] to be ~ 0.17 dex. Since the average uncertainties on the 
X-ray luminosities are of the order of 10%, the uncertainty on the 
bolometric corrections (computed using the soft and hard X-ray lu- 
minosities) is dominated by the uncertainty on L bo] and it is of the 
order of 0.20 dex. These values have to be considered as a qualita- 
tive indication on the uncertainty on L bo i and £ bo i. 



5 BLACK HOLE MASSES AND EDDINGTON RATIOS 

We estimate black hole masses from virial estimators 
jPeterson et"ai1 [2004) for 170 Type-1 AGN in our sample 
within the redshift range 0.13 ^ z ^ 3.36 (mean (z) = 1.38 
with a dispersion of .49). Of these , 96 us e the Mg II line width: 
74 sources are from Mer loni et al. 1 (120101) (with uncertainties of 



~ 0.25 dex) and 22 from lTrump et al. (2009b) (with uncertainties 
of ~ 0.4 dex). The remaining 74 are also from lTrump et al.l d2009b) 
and use the HB line width, with uncertainties of ~ 0.4 dex. We 
combine these black hole masses with bolometric luminosities to 
compute the Eddington ratio of each source: 



(2) 



L bo i L bo i 

/tEdd = x ■ 

^Edd Mbh 

Virial estimators are unavailable for the Type-2 sample; instead, 
we exploit the well-studied correlation bet ween black hole and 
host galaxy bulge mass (in particular that of lHaring & Rixll2004h 
to estimate black hole masses. We combine estimates of each host 
galaxy's total (i.e., bulge+disk) stellar mass from our SED fitting 
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Figure 5. Comparison between the values of bolometric luminosity com- 
puted by the SED-fitting code and those obtained integrating the rest-frame 
SED from 1/jm to the X-ray for the Type-1 AGN sample. Black and orange 
symbols represent the spectro-z and the photo-z sample respectively. Opens 
squares represent AGN with an upper limit in the soft X-ray band. The red 
solid line and the dashed lines represent the one-to-one correlation and 1-cr 
scatter of 0.24 dex, respectively. 



with bulge fractions from morphological assessments in order to 
determine the stellar mass of the bulge component of each host 
galaxy. This technique and the uncertainties therein are described 
below. 

We use secure morphological information for 144 Type-2 
AGN obtaine d from the updated Zurich Estimator of Structural 
Types (ZEST; IScarlata et al.ll2007h . known as ZEST+ (Carollo et 
al. 2012, in preparation). ZEST+ classifies galaxies in five morpho- 
logical types located in specific regions of the 6-dimensional space. 
Combining these measured morphologies with the results of exten- 
sive AGN host galaxy simulatio ns mapping observed m orphology 
to intrinsic bulge-to-total ratio dSimmons & Urry||2008k yields the 
following bulge fractions for each ZEST+ morphological type: 

• Elliptical: B/Tot = 1 

• S0:B/Tot = 0.5 

• bulge-dominated disk: B/Tot = 0.75 

• intermediate-bulge disk: B/Tot = 0.5 

• disk-dominated: B/Tot = 0.25 



Following lsimmons & Urrvl d2008l) and lSimmons et al.lj20TTh . un- 
certainties in B/Tot are typically ~ 0.3, with the limit that < 
B/Tot < 1. For 57 spectro-z Type-2 AGN there is significant uncer- 
tainty in the morphology (e.g., due to artifacts in the HST-ACS im- 
ages), while 47 spectro-z objects lie outside the ACS tiles; ZEST+ 
is therefore unable to determine the morphology of these sources. 
For these sources and for the photo-z sample, we employed the total 
stellar mass in the black hole mass estimate (i.e., B/Tot = 1) and 
consider the black hole mass as an upper limit. This does not affect 
our results, as we will show in Section |6l4l by treating separately 
Type-2 AGN with reliable morphological classification. 

Uncertainties in the stellar masses are mainly due to two fac- 
tors. First, each input parameter used to compute M, has an un- 



certainty: for example, the metal licity (which w e fixed to the solar 
value), the extinction law dCalzetti et al] |2000), the assumed star 
formation history, the assumed initial mass function (which is just a 
scale factor, Log M»(SalpeterIMF) ~ Log M„(ChabrierIMF)+0.23, 
we have assumed the Chabrier IMF), and different stellar popula- 
tion synthesis models. These inputs carry non-negligible uncertain- 
ties and several works in the literature have explored the impact of 
these uncertainties on the derived physical properties of galaxies 
(e.g. iMarchesini et alj|2009l ; IConrov et al]|2009l : IBolzonella etai] 
2010). The uncertainties on input parameters produce an average 
scatter of ~ 0.15 dex. 

A second source of uncerta inty is given by the d ata-set used 
and it varies for each source. IBolzonella et al.l Q2010) performed 
several tests on simulated catalogs and find a scatter of ~ 0.15 dex 
due to this effect. Combining these two sources of uncertainty, we 
consider a global uncertainty of ~ 0.2 dex on the derived values of 
M,. 

In order to properly account for these sources of uncertainty 
in the black hole mass estimates for our Type-2 AGN, we employ 
a Monte Carlo method, simulating 10 5 data points for each of our 
sources within the errors. The method also co nsiders the evolution 
of the bulge mass-black hole mass relation dMerloni et ai]|2004h 
and th e intrinsic scatter in the original relation of lHaring & Rixl 
d2004l) . We report the median mass of the Monte Carlo distribu- 
tion for each source and compute the asymmetric uncertainties 
on the mass based on the distribution. For the 481 Type-2 AGN 
with Lboi available, we further compute the Eddington ratio, addi- 
tionally accounting for uncertainties in L boI . Uncertainties in black 
hole masses and Eddington ratios are slightly asymmetric for those 
sources with asymmetric uncertainties in B/Tot, but are typically 
~ 0.5 dex for the black hole masses and ~ 0.55 dex for the Ed- 
dington ratios. The median black hole mass for the Type-2 AGN 
sample is 10 8 M Q , while the upper and lower quartiles correspond- 
ing to 75% and 25% are 2.3 x 10 s M Q and 4.8 x 10 7 M 0) respec- 
tively. The median black hole mass for the Type-1 AGN sample is 
2.7 x 10 8 M Q , while the upper and lower quartiles corresponding to 
75% and 25% are 4.3 x 10 s M B and 1.2 x 10 8 M , respectively. 

We show black hole masses and Eddington ratios for the Type- 
1 and Type-2 AGN samples in Figure [6] Note that black hole 
masses for Type 2 AGN would increase by a factor of ~ 1.7 us- 
ing a Salpeter IMF, decreasing the difference between the average 
Type-1 and Type-2 black hole masses by a similar factor. For an 
extensive discussion and comparison between Type-1 and Type-2 
AGN see Sects.|6"4l|r31and|rr6l 



6 RESULTS 

6.1 Bolometric correction vs. Bolometric luminosity for the 
Type-1 AGN sample 

We have computed the nuclear bolometric luminosities and the 
bolometric corrections for the B-band, the soft and hard X-ray 
bands considering the Type-1 AGN sample as already described 
in Sect. [4] Initially, we have computed k bol as a function of L bol 
using the sample of 373 Type-1 AGN (both spectro-z and photo-z 
samples). Subsequently, we have divided the main sample in dif- 
ferent sub-samples. We have considered all sources with both soft 
and hard X-ray detections (361 Type-1 AGN), only the spectro-z 
Type-1 AGN sample (310 Type-1 AGN) and the spectro-z sample 
but with detection in both soft and hard bands (304 Type-1 AGN). 
For each sub-sample we have estimated the relations between fcboi 
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Figure 6. Left panel: Histogram of estimated black hole masses from virial estimators for Type-1 AGN (blue filled histogram), and from stellar masses for 
Type-2 AGN (red open histogram). The two dashed lines show the median Mbh values for Type-1 AGN Mbh = 2.7 X 10 s M Q (right), and for Type-2 AGN 
Mbh ~ 10 8 M Q (left). Right panel: Histogram of estimated Eddington ratios for Type-1 AGN (blue filled histogram), and for Type-2 AGN (red open histogram). 
The two dashed lines show the median /igdd values for Type-1 AGN /Igdd = 0-12 (right), and for Type-2 AGN = 0.06 (left). Counts are normalized in 
order to better compare the samples. Black hatched histogram represents the subsample of Type-2 AGN (344 objects) for which B/Tot ratios are fixed to one 
and A/bh estimates are considered upper limits. 



and L bol in different bands. We have followed a two step procedure. 
First, we have binned the sample in order to have approximately the 
same number of sources in each bin. Then, we have computed the 
median in each bin and we have estimate d the standard devia tion in 
the median as o- mcd = 1.4826 MAD/ yfff faoaglin et all 19831) . The 
MAD term is the median of absolute deviation between data and 
the median of data (MAD = (ABS(rf - (d))), where d are the data). 
Thereby, we have fitted the median values in the six bins using a 
third-order polynomial relation, and the l<r dispersion is obtained 
using a 3.5cr clipping method. In Figure[7]the bolometric correc- 
tions as a function of the bolometric luminosity in the [0.5-2]keV, 
[2-10]keV bands and in the B-band at 0.44yum for the Type-1 AGN 
sample are presented. The spectro-z sample with detection in both 
X-ray bands is plotted with the black points, while sources with 
spectroscopic redshift and an upper limit in the soft X-ray band 
are represented with black open triangles. Orange symbols repre- 
sent the photo-z sample. Orange squares are photometric sources 
with detection in both X-ray bands, while orange open triangles 
are Type- 1 AGN that belong to the photo-z sample with un an up- 
per limit in the soft X-ray band. The red solid line represents one 
best-fit relation using the entire Type-1 AGN sample. The best-fit 
relations using different sub-samples are in close agreement. How- 
ever, in Table [2] the fcboi as a function of L bo i relations in different 
bands and for different sub-samples are reported for completeness. 
These relations approximately cover two orders of magnitudes in 
bolometric luminosities (11 < Log Lb i[L ] < 13). Sources plot- 
ted with open symbols in Fig.|7]should be considered upper limits, 
given that these L bo i are computed with an upper limit in the soft 
X-ray band. Consequently, these AGN are likely to move towards 
lower k bot and L bo i in the k ba] - L bo i plane. The effect on L bo i can- 
not be very large, unless the limit on the [0.5-2] keV fluxes is so 
low that it implies an extremely flat spectrum and therefore a total 
L x (through extrapolation) too high. Moreover, the fact that there 
is no significant difference between best-fits with and without up- 
per limits implies that neither the upper limits distribution nor the 
number of upper limits in each bin are affecting our results. The 
k bo \ - Lboi relation in the B-band seems to be nearly flat. More- 
over, the bins at the lowest bolometric luminosity decrease with 
decrease L bo i, but this could be due to several effects. First, at lower 
luminosities the statistics is poor and the decrease may be simply a 



statistical effect. Second, luminosities in the B-band might be over- 
estimated be cause of th e contr i bution of the host-galaxy emission. 
As shown in lHao et" ail <20ldl20Tl|) and Elvis et al. (2012, sub- 
mitted) from an analysis of Type-1 AGN in XMM-COSMOS the 
host-galaxy contribution in optical/near-infrared is not negligible 
and may be substantial for low luminosity AGN. Therefore, bolo- 
metric corrections at 0.44 fim for low-luminosity AGN are more 
affected by the galaxy emission making these values uncertain and 
most likely to be underestimated. 

We have used optical spectra in order to have an indepen- 
dent estimate of the possible degree of contamination by the host- 
galaxy light. We generated composite spectra for each bolometric 
luminosity bin by averaging all the available zCOSMOS spectra 
included in that bin. To create the composite, each spectrum was 
shifted to the rest-frame according to its redshift and normalized in 
a common wavelength range, always present in the observed spec- 
tral window. The composites have been then fitted using a combi- 
nation of two spectra, one representing the central active nucleus, 
the other one descri bing the host galaxy. The sets of SDSS com- 
posite spectra from iRichards et al.l d2003l) were chosen as repre- 
sentative of the quasar emiss ion, while a grid of 39 theoretical 
galaxy template spectra from iBruzual & Charlotl (120031 hereafter 
BC03), spanning a wide range in age and metallicity, were used 
to account for the stellar component. In the two lowest luminos- 
ity bins (10.2 < Log L bol < 11.8 L Q ), the zCOSMOS compos- 
ites can be fitted only if, along with an SDSS quasar composite, 
a significant host galaxy component is also included. The spec- 
troscopic host component, fitted with a bulge-dominated BC03 
template, contributes about 30% and 20% to the total luminos- 
ity at 4400A for the first bin (10.2 < Log L bol < 11.5 L Q ) and 
the second bin (11.5 < Log L DO i < 11.8 L Q ), respectively. For 
both luminosity bins, the quasar template adopted is the "dust- 
reddened" one (see Richards et al. 2003), the reddest of the com- 
posite set, suggesting that, along with host galaxy contamination, a 
fraction of our Type-1 AGN i s also experiencing a si gnificant nu- 
clear dust extinction (see also lGavignaud et al]|2006j) . In the third 
bin (11.8 < Log L bo i < 12 L Q ) the average Mb is of the order of 
-23. This value is traditionally ta ken as the threshold sep arating the 
Seyfert and quasar regimes (see Vand en Berk et al.bOOoT) . For bins 
of Log Lboi S 5 11.8 L there is no detectable host galaxy compo- 
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Figure 7. Bolometric correction as a function of the bolometric luminosity at [0.5-2]keV, [2-10]keV and in the B-band at 0.44/.im for the Type-1 AGN sample. 
The spectro-z sample with detection in both X-ray bands is represented with the black points (N=304), while sources with spectroscopic redshift and an 
upper limit in the soft X-ray band are represented with black open triangles (N=6). Orange symbols represent the photo-z sample (N=63). Orange squares 
are photometric sources with detection in both X-ray bands (N=57), while orange open triangles are Type- 1 AGN that belong to the photo-z sample with an 
upper limit in the soft X-ray band (N=6). Red bins are computed using 373 Type-1 AGN (6 bins with about 62 sources per bin). Red points are the median of 
the sources in each bin, the bars in the y axis represent the error on the median (1.4826 MAD/ VMmi). while the bars in the x axis are the width of the bin. The 
red solid line represents the best-fit relations using a third-order polynomial, while the dashed lines represent Icr dispersion after performing a 3.5cr clipping. 
Bolometric corrections have uncertainties of ~ 0.2 dex, while L(, i uncertainties are of the order of 0.17 dex, as shown by the typical error bars in the upper 
left of each panel. 



nent, and the zCOSMOS average spectrum is well fitted with an 
SDSS quasar composites alone, although again one of the reddest 
composite. 



Summarizing, the median bolometric correction of the two 
bins at bolometric luminosities less than 11.8 L e is likely to be 
a lower limit. After a proper correction of the B-band luminosity, 
bolometric corrections should increase leading to a median value 
closer to that predicted by M04 and H07. 



We did not find any relation between the bolometric correction 
drawn from a given band and the corresponding luminosity. This 
holds for both the soft and hard X-ray bands, and for the B-band 
as well. We have tested whether any relationship between k\, \ and 
luminosity exists by trying all possible permutations, but also in 
this case the data distribution is flat (see Vasudevan & Fabian l2007l 
and LI 1 for similar results). 



6.2 Bolometric correction vs. Bolometric luminosity for the 
Type-2 AGN sample 

The same analysis presented in the previous Section has been ap- 
plied to the Type-2 AGN sample. From the main sample of 547 
Type-2 AGN, 488 AGN (~ 89%) have an estimate of the bolomet- 
ric luminosity, and therefore an estimate of the bolometric correc- 
tion is available from the SED-fitting code. For obscured AGN the 
optical emission is mostly dominated by the host-galaxy, hence we 
cannot estimate the nuclear luminosity in the B-band at 0.44yum. 
The 488 objects have been divided in subsamples as already done 
for Type-1 AGN. The sample is composed by 180 spectroscopic 
Type-2 AGN with both X-ray detections, 68 objects with spectro- 
z and [0.5-2]keV upper limits, 161 photo-z Type-2 AGN with both 
X-ray detections, and 79 photo-z Type-2 AGN with [0.5-2]keV up- 
per limits. In Figure [8] the bolometric corrections as a function of 
the bolometric luminosity in the [0.5-2] and [2-10]keV bands for 
the Type-2 AGN sample are presented. These relations cover more 
than two orders of magnitudes (10 < Log L bo i[L G ] < 12.5).Also 
for the Type-2 AGN sample, the best-fit relations using different 
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Table 2. Bolometric correction relations in different bands for the X-ray selected Type-1 and Type-2 AGN samples. 

Sample N a\ ai fo" <r N b Band 

x x 2 x 3 3.5o" clipping 



Type-1 







0.239 


0.059 


-0.009 


1.436 


0.26 


1 


[0.5 - 2]keV 


vnppfrn -U Prints 


37 3 


288 


111 


-0 007 


1 308 


26 


[ 








-0.011 


-0.050 


0.065 


0.769 


0.23 


6 


B-band (0.44 /im) 






0.248 


0.061 


-0.041 


1.431 


0.26 


1 


[0.5 - 2]keV 


Spectro + Photo 


361 


0.310 


0.114 


-0.020 


1.296 


0.25 


2 


[2 - 10]keV 


[0.5 - 2]keV detected 




-0.026 


-0.037 


0.075 


0.760 


0.22 


5 


B-band (0.44 /im) 






0.250 


0.044 


-0.023 


1.455 


0.26 


1 


[0.5 - 2]keV 


Spectro 


310 


0.261 


0.108 


0.009 


1.331 


0.25 


1 


[2 - 10]keV 






0.041 


-0.065 


0.028 


0.763 


0.22 





B-band (0.44 /im) 






0.219 


0.068 


0.007 


1.444 


0.25 


2 


[0.5 - 2]keV 


Spectro 


304 


0.245 


0.123 


0.024 


1.321 


0.25 


1 


[2 - 10]keV 


[0.5 - 2]keV detected 




0.023 


-0.051 


0.043 


0.754 


0.22 





B-band (0.44 /im) 


Type-2 






0.217 


0.009 


-0.010 


1.399 


0.27 


2 


[0.5 - 2]keV 


Spectro + Photo 


488 


0.230 


0.050 


0.001 


1.256 


0.25 


2 


[2 - 10]keV 


Spectro + Photo 




0.208 


-0.059 


-0.038 


1.455 


0.28 





[0.5 - 2]keV 


[0.5 - 2]keV detected 


341 


0.217 


-0.022 


-0.027 


1.289 


0.26 





[2 - 10]keV 






0.293 


0.0652 


0.0029 


1.470 


0.25 


1 


[0.5 - 2]keV 


Spectro 


248 


0.386 


0.071 


-0.010 


1.395 


0.23 


1 


[2 - 10]keV 


Spectro 




0.275 


0.104 


0.017 


1.459 


0.25 





[0.5 - 2]keV 


[0.5 - 2]keV detected 


180 


0.411 


0.086 


-0.010 


1.395 


0.23 





[2 - 10]keV 



[a] y = ai x + 02 x ~ + a 3 x 3 + b, where x = Log L - 12 (L is the bolometric luminosity in units of L Q ) and y = Log [L/iBand]- 

[b] Number of objects removed from the cr clipping method. 



sub-samples are not significantly different. In Table [2]the £ bol as a 
function of L bol relations in the X-ray bands and for different sub- 
samples are reported. 

6.3 Bolometric correction vs. Bolometric luminosity: 

comparison of the results for Type-1 and Type-2 AGN 

Figure|9]shows the bolometric correction as a function of the bolo- 
metric luminosity in the [0.5-2]keV, [2-10]keV bands and in the 
B-band with lcr dispersions after performing a 3.5cr clipping, for 
the Type-1 AGN sample (N=373) and for the Type-2 AGN sam- 
ple with AGN best-fit (N=488), respectively. As a comparison, the 
predicted curves obtained by H07 and M04 in the soft, hard X-ray 
bands and in the B-band with lcr dispersion are also reported. De- 
spite the large scatter, the trend of increasing bolometric correction 
at increasing bolometric luminosity is confirmed. Moreover, it is 



evident that the k bol - L boi relations by M04 and H07 are higher 
than those observed in our (X-ray selected) samples for both Type- 
1 and Type-2 AGN. The different normalization of the k bo i - L bo \ 
relation between H07 and M04 has to be ascribed to a different def- 
inition of bolometric luminosity. H07 include infrared wavelengths 
as they are interested in an empirical bolometric correction, while 
fcboi values in M04 are estimated neglecting the optical-UV-X-ray 
luminosity reprocessed by the dust and therefore are representative 
of the AGN accretion power. In the present analysis we have also 
considered only the accretion powered luminosity and thus our val- 
ues should be compared with the M04 curves. The normalization 
shift between our values for the A: bo i - L bo i relationships and the 
M04 curves could be, at least partly, due to a selection effect. Since 
our sample is hard X-ray selected, it is biased towards X-ray bright 
objects with lower bolometric corrections. The SED adopted by 
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Figure 8. Bolometric correction as a function of the bolometric luminosity in the [0.5-2]keV and [2-10]keV bands for the Type-2 AGN sample with AGN 
best-fit. Symbols key are as in Fig. [7] The sample used to compute the bins is composed as follow: 180 spectroscopic Type-2 AGN with both X-ray detections 
(black points), 68 objects with spectro-z and [0.5-2]keV upper limits (black open triangles), 161 photo-z Type-2 AGN with both X-ray detections (orange 
squares), 79 photo-z Type-2 AGN with [0.5-2]keV upper limits (orange open triangles). Red bins are computed using the 488 Type-2 AGN sample (6 bins 
with about 80 sources per bin). 




Log L bol /L 

Figure 9. Bolometric correction as a function of the bolometric luminosity in the soft, hard and B-bands for the Type-1 AGN sample (N=373, black solid 
line), soft and hard bands for the Type-2 AGN sample with AGN best-fit (N=488, red solid line). The lcr dispersion after performing a 3.5<x clipping is also 
reported with the dashed lines. The green and blue lines represent the bolometric correction in the hard band and lcr dispersion obtained by H07 and M04, 
respectively. 
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M04 is typical of luminous optically selected quasars, therefore it 
is biased towards higher k bol values. 

The present results may have interesting consequences. In fact 
all accretion models, that also include mergers, fail in reproduc- 
ing the high mass e nd of the local black hole mass function (see 
Shan kar et al] |2009). However, as suggested by our data, assum- 
ing a lower k bo \ than the one inferred by M04 or H0 7, can ease the 
tensio n between models and data (see discussion in Shan karet al.l 
l20lTh . The range of validity of these curves are limited to slightly 
more than two orders of magnitudes for both Type-1 and Type-2 
AGN, where L bol ranges from 10 10 L o to 10 12 L Q for Type-2 AGN, 
and from 10"L o to 10 13 L o for Type-1 AGN. In the overlapping 
range of bolometric luminosity (10 n ~ 12 L Q )> there is no significant 
difference between the bolometric corrections for Type- 1 and Type- 
2 AGN. It is also interesting and noteworthy that the k bol - L bol 
relations for Type-2 AGN seem to be the natural extension of the 
Type-1 relations at lower luminosities. Even if we can explore a 
limited range of bolometric luminosities in both AGN samples this 
relationship can be applied for all AGN across nearly four decades 
in luminosity. 

As a final comment, we want to discuss a comparison be- 
tween the results presented in this paper and the results on the 
a ox parameter in L10. The fact that a ox and k bo \ show a simi- 
lar behaviour is a natural consequence of the tight correlation be- 
tween these two parameters, which has been discussed in depth 
in L10. Indeed, a ox is almost independent on the X-ray luminos- 
ity at 2 keV, while there is a strong trend with the optical lumi- 



nosity at 2500A (see alsolTananbaum et aQ 1979 ; Zamorani et al.l 
198 1[ IVignali et alj [20031; Isteffen et alj l2006r . ljust et all 12007k 



Young et alj|2009l ; iMarchese et alj|20i2h . which is a tracer of the 



bolometric luminosity (more than 70% of L bo \ comes from the 
optical-UV). It is not surprising then to find a correlation between 
k bol (evaluated in the soft and the hard X-ray bands) and L bo i, while 
no correlation is seen with the X-ray luminosity. In conclusion, this 
analysis suggests that the fundamental underlying correlation arises 
between k hol and L bo[ . 

6.4 Hard X-ray bolometric correction vs. Eddington ratio 

Several works in the literature found a trend between the hard 
X-ray fc bo i and /i Edd , although wi th the presence of a large scat- 
ter (e .g., IVa"sudevan & Fabianll2007l ; V09a; V09b. lVasudevan et al.l 
l20ld . L10), or betwee n i Edd and th e intrinsic bolometric AGN lu- 
minosity (e.g., lTrumpetal.ll2009al l 201 lh . The scatter is not re- 
duced even considering AGN with simultaneous optical, UV and 
X-ray data retrieved from the XMM-Newton EPIC-pn and Optical 
Moni tor (OM) archives (see Fig. 11 in L10. IVasudevan &~ Fabian 
2009). In Figure [I0]the hard X-ray k bol as a function of the ,l Edd for 
Type-1 AGN is presented. We have computed the ordinary least- 
squares (OLS) bisector for the k bo \ - /l Edd relation considering the 
170 Type-1 AGN with M BH mass estimates from broad lines (there 
are only two objects with an upper limit in the soft band). The best- 
fit parameters for the k bo i - /l Edd relation using OLS(Y|X) (i.e. treat- 
ing /i Edd as the independent variable), OLS(X|Y) (i.e. treating the 
hard X-ray k bol as the independent variable) and the OLS bisector 
are reported in Table [3] We find that the slope of the bisector rela- 
tion considering the 170 Type-1 AGN sample is statistically consis- 
tent with the slope found considering the subsample of 150 Type-1 
AGN in L 10. 

The same analysis has been performed using 488 Type-2 AGN 
for which bolometric luminosities and stellar masses are available 
from the SED-fitting, and for the subsample of 144 Type-2 AGN 
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Figure 10. Hard X-ray bolometric correction versus Eddington ratio for the 
170 Type-1 AGN (black points) with BH mass estimate from broad lines. 
The short-dashed black line shows the best-fit relation and lcr dispersion 
that we found using the OLS bisector algorithm. Typical en'or bars in the 
upper left of the panel are showed for both /t Edd (0.43 dex) and k bo \ (0.2 
dex). 
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Figure 11. Hard X-ray bolometric correction versus Eddington ratio for the 
488 Type-2 AGN with bolometric luminosities and stellar masses available 
from the SED-fitting. Symbols key are as in Fig. [7] The short-dashed black 
line shows the OLS bisector and lcr dispersion for the 488 Type-2 AGN 
with L bo \ and M„ available. The long-dashed black line shows the OLS 
bisector and lcr dispersion for the 144 Type-2 AGN with L bo i, M„ and mor- 
phology classification available. Typical error bars in the upper left of the 
panel are showed for both /t Edd (0.55 dex) and k bo [ (0.2 dex). 
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Figure 12. Bolometric luminosities as a function of black hole masses for 
both Type-1 (blue symbols) and Type-2 (red symbols) AGN samples. Diag- 
onal lines represent the trend between L(, | and Mbh at different fractions 
of the Eddington luminosity: indd = 1,0.1 and 0.01 (solid, long-dashed 
and short dashed lines, respectively). Black open squares mark 144 Type-2 
AGN with Lij i, M„ and morphology classification available. Typical error 
bars in the lower right of the panel are showed for both Mbh (0.4 dex for 
Type-1 AGN in blue and 0.5 dec for Type-2 AGN in red) and k^ a \ (0.2 dex 
for both Type-1 and Type-2 AGN). 



Table 4. Hard X-ray bolometric correction as a function of the Eddington 
ratio for the X-ray selected Type-2 AGN sample. 



(m ± dvrif 



(q ± dqY cr [3.5<r clipping] N Relations 



488 Type-2 C 


0.621 ± 0.025 


1.947 ±0.035 


0.34 


2 


OLS bisector 


0.283 ± 0.036 


1.536 ±0.052 


0.26 


4 


OLS(Y|X) 


1.105 ±0.084 


2.536 ±0.104 


0.57 





OLS(X|Y) 


144 Type-2 d 


0.649 ± 0.054 


2.090 ± 0.084 


0.32 





OLS bisector 


0.375 ± 0.097 


1.702 ±0.151 


0.23 


3 


OLS(Y|X) 


1.015 ±0.095 


2.608 ±0.137 


0.44 


1 


OLS(X|Y) 



[a] y = (m ± dm) x + (q ± dq), where x = Log /Igdd and 
v = Log [L/L[2-10]keV]- 

[b] Number of objects removed from the cr clipping method. 

[c] Type-2 AGN with L^oi and M» available. 

[d] Type-2 AGN with L(, ], M» and morphology classification available. 



Table 3. Hard X-ray bolometric correction as a function of the Eddington 
ratio for the X-ray selected 170 Type-1 AGN sample with Mbh available. 



(m ± dmf 


(q ± dqf 


cr [3.5cr clipping] 


Relations 


0.752 ± 0.035 


2.134 ±0.039 


0.32 


OLS bisector 


0.392 ± 0.055 


1.799 ±0.063 


0.27 


OLS(Y|X) 


1.303 ±0.149 


2.647 ±0.128 


0.49 


OLS(X|Y) 



[a] y = (m ± dm) x + {q± dq), where x = Log /lEdd and 
y = Log [L/L[ 2 -io]kev]- 

No object has been removed from the cr clipping method. 

with reliable morphology classification. The data are shown in Fig- 
ureQj] The trend of increasing Eddington ratios at increasing bolo- 
metric corrections, as found for Type-1 AGN, is confirmed also for 
Type-2 AGN. The slope of the bisector relation considering the 170 
Type-1 AGN sample is marginally consistent at the 3<x level with 
the slope found considering the 488 Type-2 AGN sample, while it 
is fully consistent with the slope found considering the 144 Type-2 
AGN sample with reliable morphological classifications. Also the 
normalizations of the kb a \ - /Irad relations for the Type 1 and Type 2 
AGN are in good agreement with each other. The best-fit parame- 
ters for the fcboi - /iEdd relation using OLS(Y|X), OLS(X|Y) and the 
OLS bisector are reported in Table [5] 

The hard X-ray radiation of AGN with relatively high Edding- 
ton ratio is commonly thought to be produced from a disk corona 
as a result of Comptonization of soft photons arising from the ac- 



cretion disk (e.g lHaardt & Maraschi|[l99ll 1 19931 : iKawaguchi et al.l 
l200ll; lcaol[2009). If the bolometric luminosity is the result of ac- 
cretion disk and corona emission, the fraction of X-ray luminosity 
over the total luminosity represents the strength of the corona rel- 
ative to the accretion disk. Therefore, the correlation between the 
hard X-ray bolometric correction and the Eddington ratio for Type- 

1 and Type-2 AGN indicates that the corona relative to the disk 
becomes weaker as the Eddington-scaled accretion rate increases. 

In Figure [T2]the bolometric luminosities are plotted as a func- 
tion of black hole masses for both Type-1 and Type-2 AGN sam- 
ples. In the following we will focus on the sources which are un- 
doubtedly dominated by AGN activity (i.e. we have removed 7 
sources from the Type-2 AGN sample with Log L[2-iojkeV < 42 
erg s~'). Diagonal lines represent the trend between L bol and M Bli 
at different fractions of the Eddington luminosity (A EM = 1,0.1 
and 0.01). It is worth noting that black hole masses for Type-1 and 
Type-2 AGN are derived in completely different ways (i.e., virial 
estimators versus scaling relations), and the difference between 
black hole masses (and Eddington ratios) for Type-1 and Type- 

2 AGN would decrease if a Salpeter IMF were used to compute 
stellar masses for Type-2 AGN. Stellar masses computed with the 
Salpeter IMF would increase by a factor ~1.7. As a consequence, 
black hole mass and Eddington ratio estimates would increase by a 
similar factor. There is a continuity between bolometric luminosi- 
ties as a function of Mbh for Type-1 and Type-2 AGN, where few 
sources have L bo] lower than 0.01 L Ed d- To determine whether Ed- 
dington ratios are affected by any significant evolution, we have 
studied a possible dependence of A^d on redshift, Lboi, X-ray lu- 
minosities, M B h an d column densities. We find no correlation be- 
tween Eddington ratios with both X-ray luminosities and column 
densities. The other correlations are discussed in the following. 
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6.5 Luminosity-redshift dependence of Eddington ratio 

We have explored the possibility of a dependency of A EAi with red- 
shift and bolometric luminosity by binning both Type-1 and Type-2 
AGN samples in z and L bol . The samples are divided in two red- 
shift bins and three Lboi bins. The redshift bins are z < 1.2 and 
1.2 < z < 2.3, while the luminosity cuts for each sample are cho- 
sen in order to have approximately the same number of objects in 
each bin. The redshift bins have been defined in order to sample the 
observed evoluti on of the hard X -ray luminosity function of AGN 
determined bv lAird et al.ld2010l see their Fig. 9 and the discussion 
below). There are 60 Type-1 AGN and 317 Type-2 AGN at z < 1.2, 
while 109 Type-1 AGN and 135 Type-2 AGN are at 1.2 < z < 2.3. 
At z > 2.3 the number of AGN is not large enough to be statistically 
significant. The black histograms in Figs. 1 131 and 1 141 show the ob- 
served A EM distributions for Type-1 and Type-2 AGN, respectively. 

The observed Aem distributions may be biased by selection ef- 
fects related to the depth of the X-ray data, and the fall-off below 
the peak at low A EiA can be partly due to incompleteness of the 
X-ray selection. We have therefore quantified the impact of this in- 
completeness on our A EM distribution f or each bin by employing the 
standard V max method, introduced by Schmidt] dl968r) . The quantity 
Vmax represents the maximum volume where an object would still 
be detectable in our survey given its X-ray luminosity, redshift and 
column density and is described by 



43-45.03 



45.03-45.48 45.48-46.5 



^max ( 



r 



CXL x (i),z(i),N H (Ml + zf^dz, (3) 
dz 



where z min is the lower boundary of the redshift bin, z max is the 
minimum between the upper boundary of the redshift bin and the 
redshift where the i th object would not be detectable anymore in 
the survey. The parameter CI is the solid angle covered by XMM- 
Newton at the flux level f(Lx(i), z), and dV/dz is the comoving vol- 
ume. The term (1 +z) k describes the AGN evolution, which we have 
chosen to represent as a pure density evolutioifl At zeroth-order, 
we have neglected any Nh dependence and, therefore, we have es- 
timated the flux at each redshift considering a simple k-correction 



/(ix(0,z) 



L x (i)(l + z) ([ - a) 



Andi 



(4) 



where a = 0.7 dLa Franca et al.ll2005h and L\ is the de-absorbed 
X-ray luminosity in the [2 - 10] keV rest-frame band. To obtain 
the area at each X-ray flu x we have employed t he [2 - 10] keV 
sky coverage computed by Cappelluti et al. (2009, see their Fig. 5), 
and we have finally integrated over the comoving volume. The k 
value in the evolutionary term has been chosen to match the ob- 
served evolu tion of the hard X -ray luminosity function of AGN de- 
termined bv lAird et al.l d2O10l see their Fig. 9). In the low redshift 
bin we have considered k = (2.2, 2.4, 4.3, 7.8) for Log Lp-ioikev = 
(42,43,43.5,44) erg s" 1 , while for Log L [2 -io]keV > 44 erg s _1 we 
have used 7.8^. We have adopted k = consistently with the lumi- 
nosity function at z > 1 .2 by Aird et al. for the range of luminosities 
covered by our sample. For each AGN we have adopted the appro- 
priate k value interpolating the numbers above. The total volume 
(V) has been estimated through the same procedure, where now 
the area is the total area covered by XMM-Newton (i.e. 2.13 deg 2 ). 
The A EM distributions weighted by the ratio V/ V raax for each object 
are plotted with the red dashed histograms in Figs.[l3]and[T4] As 



The results would be similar if we adopted a different evolutionary law. 




-3 -2 -1 1-3 -2 -1 1-3 -2-10 1 



Loe A,. 



Log A E 



The choice of this value is not critical 
Log L [2 -io]keV > 44 erg s" 1 have V/V max =s 1. 



because all objects at 



Figure 13. Distributions of Eddington ratios in bins of luminosity and red- 
shift for the Type-1 AGN sample. The panels are divided between z. < 1.2 
(top three panels) and 1,2 < z < 2.3 (lower three panels), and bolometric 
luminosities increases from left to right (Log Lboi [erg s~'] intervals are re- 
ported on top of each panel). The black histograms show the observed A E ^ 
distributions, while the red dashed histograms represent the distributions 
corrected for incompleteness, according to the V*/V* ma x method described 
in § 16.51 The dashed red lines are the median values corresponding to the 
completeness-corrected distributions, while the solid black lines in the first 
and in the second bins are plotted at the A E( n value corresponding to the 
median in the highest luminosity bin. There are ~ 20 and ~ 37 objects in 
each bin at z < 1.2 and 1.2 < z < 2.3, respectively. The numbers of AGN in 
the completeness-corrected A E ^ distributions are 98, 22, and 21 at z < 1.2, 
while there are 88, 51, and 38 at 1.2 < z < 2.3, from left to right. 



shown in the figures, the completeness correction does not change 
significantly the histograms. There are several interesting points to 
note. 

First, the distribution of Log A EM is nearly Gaussian, espe- 
cially at high redshift and high luminosity, with a dispersion of the 
order of ~ 0.35 dex for the Type-1 AGN sample and ~ 0.5 dex for 
the Type-2 AGN sample. As expected, the low redshift/luminosity 
bins are more sensitive to incompleteness. 

Second, it is evident that the population of AGN that we are 
studying is dominated by sub-Eddington accreti on rate objects. 
This re sult in contrast with the findings obtained bv lKollmeier et al] 
(2006), where the AGN population is dominated by near-Eddington 
accretors. The difference, consistently with the trend that we see 
in our data (see below), might be due to the fact that the bulk of 
the AGN population studied by Kollmeier and collaborators have 
higher L bo i, typically in the range Lboi ~ 10 45 ~ 47 erg s~'. 

Third, the Eddington ratio increases with luminosity for both 
Type-1 and Type-2 AGN. In Fig. [T5] the median A E a is plotted 
against the median L bol for both AGN populations. Different sym- 
bols for low redshift (filled circles) and high redshift (open squares) 
are introduced. There is no clear evolution of A EiA with redshift 
in both redshift bins. At a given L bo] , Type-2 AGN seem to have 
higher A E a than Type-1 AGN at low redshift, while at high red- 
shift the difference is not significant. A summary of the average 
Log A Eii values and relative dispersions for the Type- 1 and Type-2 
AGN samples are given in Tables [5] and [6] respectively. 

Summarizing, we have found that Eddington ratio evolves 
with bolometric luminosity for both Type-1 and Type-2 AGN, 
while it does not show a clear evolution in redshift if we bin in 
i-boi- Type-1 AGN have median Eddington ratios ranging, on av- 
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Figure 14. Distributions of Eddington ratios in bins of luminosity and red- 
shift for the Type-2 AGN sample. Description as in Fig. 1131 There are ~ 106 
and ~ 45 objects in each bin at z < 1.2 and 1.2 < z < 2.3, respectively. The 
numbers of AGN in the completeness-corrected /l Edd distributions are 1317, 
234, and 140 at z < 1.2, while there are 259, 77, and 62 at 1.2 < z < 2.3 
from left to right. 



Table 5. Averaj 


;e Log /l Edd values in L[,ol - z bins for Type-1 AGN. 




L °§ L bol 

[erg s~'] 


(Log L bol > b 
[erg s~'] 




N d 


(Log /l Edd > 
mean 


cr 


(Log /l Edd > 
median 


z< 1.2 


43.00-45.03 


44.80 


20 


98 


-1.47 


0.38 


-1.66 


45.03-45.48 


45.30 


20 


22 


-1.15 


0.33 


-1.18 


45.48-46.50 


45.63 


20 


21 


-0.78 


0.47 


-0.83 


1.2<z« 2.3 


44.6-45.53 


45.35 


37 


88 


-1.01 


0.37 


-1.04 


45.53^15.95 


45.74 


36 


51 


-0.83 


0.31 


-0.84 


45.95-48.00 


46.22 


36 


38 


-0.56 


0.35 


-0.58 



[a] Bolometric luminosity intervals. 

[b] Median Bolometric luminosity. 

[c] Number of observed objects in each bin. 

[d] Number of predicted objects in each bin employing the V/V mat method. 




44 44.5 45 45.5 46 

Lo g L Bol [ er S s_1 ] 

Figure 15. Median Eddington ratios as a function of the median Lb„i for 
Type-1 (blue symbols) and Type-2 AGN (red symbols). Filled circles and 
open squares represent the median /l Edd for z < 1.2 and 1.2 < z < 2.3, 
respectively. Solid lines connect low redshift bins, while dashed ones con- 
nect high redshift bins. Error bars on the median are estimated considering 
the 16th and 84th percentile divided by the square root of the number of 
observed AGN. 



erage, from Log A Eid ~ -1.6 to Log /l Edd ~ -0.6 across the lu- 
minosity scale (with a dispersion of ~ 0.35 dex), while the corre- 
sponding values for Type-2 AGN range from Log A EM 1.8 to 

Log /l Edd 0.7 (with a dispersion of ~ 0.5 dex). 



Table 6. Average Log /l Edd values in L(, i - z bins for Type-2 AGN. 




[erg s~'] 


(Log L bo i> b 
[erg s~'] 


N ob S 


N d 

coir 


(Log i Edd ) 
mean 


cr 


(Log i Edd > 
median 


z < 1.2 


43.00-44.35 


44.08 


106 


1317 


-1.72 


0.59 


-1.83 


44.35^14.81 


44.62 


106 


234 


-1.41 


0.50 


-1.49 


44.81^16.80 


45.08 


105 


140 


-0.97 


0.57 


-1.08 


1.2 < z s? 2.3 


44.30-45.10 


44.88 


44 


259 


-1.25 


0.55 


-1.24 


45.10-45.52 


45.35 


45 


77 


-1.01 


0.43 


-0.97 


45.52-48.00 


46.22 


46 


62 


-0.71 


0.46 


-0.74 



[a] Bolometric luminosity intervals. 

[b] Median Bolometric luminosity. 

[c] Number of observed objects in each bin. 

[d] Number of predicted objects in each bin employing the V/V max method. 

6.6 Black hole mass-redshift dependence of Eddington ratio 

The distribution of Eddington ratio as a function of black hole mass 
and redshift delivers more significant constraints on the physical 
distribution of the fueling rates. The observed and V/Vmax cor- 
rected Aem distributions at a given M^h and redshift are plotted in 
Figs. [T6] and [T7] for the Type-1 and Type-2 AGN sample, respec- 
tively. For both samples, each bin contains almost the same num- 
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Figure 16. Distributions of Eddington ratios in bins of black hole mass 
and redshift for the Type-1 AGN sample. The panels are divided between 
z < 1.2 (top three panels) and 1.2 < z < 2.3 (lower three panels), and black 
hole masses increase from left to right (Log Mbh [M q ] intervals are re- 
ported on top of each panel). The black histograms show the observed /l^dd 
distributions, while the red dashed histograms represent the distributions 
corrected for incompleteness. The dashed red lines are the median values 
corresponding to the completeness-corrected distributions, while the solid 
black lines in the first and in the second bins are plotted at the A E m value 
corresponding to the median in the highest luminosity bin. There are ~ 20 
and ~ 37 objects in each bin at z < 1.2 and 1.2 < z < 2.3, respectively. The 
numbers of AGN in the completeness-corrected A E ^i distributions are 96, 
26, and 19 at z < 1.2, while there are 79, 52, and 45 at 1.2 < z < 2.3 from 
left to right. 



ber of sources. From Fig. |T6]it is evident that, for Type-1 AGN, the 
completeness limit affects the A EM distribution only at low black 
hole masses in the low redshift bin. The situation for Type-2 AGN 
seems to be more complicated, where the low redshift bins are more 
affe cted by incomplete ness in all M BH intervals (see Fig.ll7t, 

lAird et all febl2h recently claimed that, for any given stellar 
mass, the Eddington ratio distribution of X-ray selected obscured 
AGN at Log L x > 42 and z < 1.0 is well described by a power law. 
They reach this conclusion by running a variety of Monte Carlo 
simulations to correct their sample for a number of incomplete- 
nesses. Given the correlation between M, and Mbh, the Aird et al. 
power-law distribution should be seen as a power-law also in bins 
of M BH . As shown in Figure[T6]and[T7] we do not see any evidence 
for such a distribution in most of our redshift and Mbh bins. This is 
particularly clear for the Type-1 AGN (not included in Aird's anal- 
ysis). The only sub-sample of Type 1 AGN, where a distribution 
continuously increasing toward low /tedd is seen, is the sample in the 
low redshift and low Mbh bin. All the other sub-samples of Type-1 
AGN show A EM distrib utions more consistent wi th Gaussians (for 
a similar result see also lSteinhardt & Elvisl l2010). The situation is 
less clear-cut for the Type-2 AGN, where there are at least two bins 
at low redshift where the completeness-corrected distributions may 
suggest the presence of an underlying distribution increasing to- 
ward low A EM . However, no clear evidence for such a power law 
is present in our data in the higher redshift bin for Type 2 AGN. 
This result is consistent with the findings of Shankar et al. (2011) 
where it is shown that the A EM distribution at high redshift has to be 
Gaussian in order to match the observed luminosity function, while 
at low redshift the power-law distribution is preferred. 

The Eddington ratio as a function of Mbh is plotted in Fig. 1181 




-3 -2 -1 1-3 -2 -1 1-3 -2-10 1 
Log A Edd Log A Edd Log A Edd 

Figure 17. Distributions of Eddington ratios in bins of M B h and redshift 
for the Type-2 AGN sample. Description as in Fig. 1161 There are ~ 106 
and ~ 55 objects in each bin at z < 1.2 and 1.2 < z < 2.3, respectively. The 
numbers of AGN in the completeness-corrected A E a distributions are 1046, 
384, and 262 at z < 1.2, while there are 201, 132, and 65 at 1.2 < z < 2.3 
from left to right. 



for the Type-1 and Type-2 AGN samples. The two AGN sam- 
ples show higher A E a at higher r e dshift at any given Mbh- In 
Shank ar et alT ll2004» and [shankar (2009) it was shown that an 
increasing A EM with redshift may yield better results with the 
low mass end of the local BH mass function. Shankar et al. 
(2011) also showed that an increasing A E a with redshift yields 
very good agreement with the high duty cyc l es inf erred from 
X-ray studies at z - (e.g., iGoulding etaD 2010h. and w ith 
[Om] lines (e.g.. iKauffmann et alj|2003l and iBest et al]|2005l) . A 
trend of increasing /i Ed d wit h redshift has also been also found by 
iNetzer & Trakhtenbrod d2007l) using a sample of 9818 SDSS Type 
1 AGN at z < 0.75. A comparison between their sample and ours is 
difficult given that we are sampling a different redshift range (only 
18 Type-1 AGN have z < 0.75 in our sample). However, we are 
in agreement with the result by Netzer and collaborators extending 
the analysis to higher redshifts and using a sizable X-ray selected 
Type-1 AGN sample. 



7 SUMMARY AND CONCLUSIONS 

A homogeneous analysis of the bolometric output and Eddington 
ratio of 929 AGN at different X-ray absorption levels is presented. 
Several aspects of the present analysis have been improved with 
respect to L10 and Lll. In particular, the far-infrared emission is 
now better constrained thanks to the inclusion of Herschel data at 
100 fim and 160 fim in the SED-fitting code for Type-2 AGN. Our 
main sample is further extended at fainter magnitudes with the ad- 
dition of a sizable number of objects with photometric redshift, in 
order to take bias and selection effects under control. The photo- 
metric redshift catalog is the latest release provided by S 1 1 . More- 
over, we have increased the coverage in the near-infrared includ- 
ing the H band photometry. Black hole masses for Type-1 AGN 
are available for 170 sources computed from virial estimators us- 
ing different lines width (Mgn and H/?). Black hole masses and 
Eddington ratios for Ty pe-2 AGN are estim ated for 481 objects 
through scaling relations (Haring & Rix 2004) using a Monte Carlo 
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Log M BH [M ] 

Figure 18. Median Eddington ratios as a function of the median Mbh for 
Type-1 (blue symbols) and Type-2 AGN (red symbols). Filled circles and 
open squares represent the median /Iridd for z < 1.2 and 1.2 < z < 2.3, re- 
spectively. Solid lines connect low redshift bins, while dashed ones connect 
high redshift bins. Error bars on the median are estimated considering the 
16th and 84th percentile divided by the square root of the observed number 
of AGN. 



method in order to account for uncertainties in M„, L bol , as well as 
the intrinsic scatter in the M BH - M, relation. We have analysed 
the dependence of k bo \ on L bo \ in the B-band at 0.44 //m, in the soft 
and hard X-ray bands and we have compared our results with the 
predicted curves by M04 and H07. Eddington ratios are studied as 
a function of hard X-ray luminosities, Mbh, and redshift for 
both Type- 1 and Type-2 AGN samples taking into account incom- 
pleteness effects. 

Our main results are the following: 

(i) There is a trend for higher bolometric corrections at higher 
bolometric luminosities in the [0.5-2]keV and [2-10]keV bands for 
both Type-1 and Type-2 AGN samples. The kboi — £*oi relations by 
M04 and H07 predict higher k hol than what is observed in our X- 
ray selected AGN samples (both Type-1 and Type-2 AGN). The 
range of validity of these curves are limited to about two orders of 
magnitudes for both Type-1 (11 < Log L bo ][L Q ] < 13) and Type- 
2 AGN (10 < Log L bo i[L ] ^ 12). In the overlapping luminosity 
range (Lb i = 10"~' 2 L G ) there is no significant difference between 
k bol for Type- 1 and Type-2 AGN, and moreover, the k bol - L bol rela- 
tion for Type-2 AGN seem to be the natural extension of the Type- 1 
relation at lower luminosities. 

(ii) The k bol - L bol relation in the B-band is in agreement with the 
H07 and M04 relations for L bo i > 6x 10 n L o . The lowest bolometric 
luminosity bins do not follow the predicted trend. A fraction of our 
Type-1 AGN is experiencing a significant nuclear dust extinction, 
along with host galaxy contamination (Elvis et al 2012, submitted), 
and these two effects are sufficient to explain the observed differ- 
ence. 

(iii) We confirm the trend to have higher k bol for higher Edding- 
ton ratios for both Type-1 and Type-2 AGN. The same trend has 
been observed with the bolometric luminosity. This indicates that 



the emission from the X-ray corona becomes weaker relative to the 
disc as the Eddington-scaled accretion rate increases. 

(iv) The population of AGN is dominated by sub-Eddington ac- 
cretion rate objects at a given L bol . 

(v) The distribution of Log is nearly Gaussian especially at 
high redshift and at high L bo i/M BH , with a dispersion of the order 
of ~ 0.35 dex for the Type-1 AGN sample and ~ 0.5 dex for the 
Type-2 AGN sample. As expected, the low redshift/luminosity bins 
are more affected by incompleteness. 

(vi) Eddington ratio increases with bolometric luminosity for 
both Type-1 and Type-2 AGN. 

(vii) Eddington ratios show an evolution in redshift if we bin in 
M BH for both AGN Types. If we instead bin in bolometric luminos- 
ity /lEdd do not show any clear evolution in redshift for both AGN 
Types. 

We want to emphasize that the k bo] - L bol relations derived in the 
present work are calibrated for the first time against a sizable AGN 
population, and therefore rely on observed redshifts, X-ray lumi- 
nosities and column density distributions. The application of these 
empirical relations offers the opportunity of future developments 
along several lines of investigation. For example, they could pro- 
vide important hints for the computation of the black hole mass 
density and AGN bolometric luminosity function. As a final com- 
ment, this analysis suggests that the fundamental physical corre- 
lation of k bol is with bolometric luminosity and Eddington ratio, 
rather than with single band luminosities. 



ACKNOWLEDGMENTS 

In Italy, the XMM-COSMOS project is supported by ASI-INAF 
grants 1/009/10/0, 1/088/06 and ASI/COFIS/WP3110 1/026/07/0. 
In Germany the XMM-Newton project is supported by the Bun- 
desministerium fur Wirtshaft und Techologie/Deutsches Zentrum 
fur Luft und Raumfahrt and the Max-Planck society. M. Salvato ac- 
knowledges support by the German Deutsche Forschungsgemein- 
schaft, DFG Leibniz Prize (FKZ HA 1850/28-1). F. Shankar ac- 
knowledges support from a Marie Curie grant. The author acknowl- 
edges the anonymous reviewer who provided many useful sugges- 
tions for improving the paper. The entire COSMOS collaboration 
is gratefully acknowledged. 



REFERENCES 

Aird J., Coil A. L., Moustakas J., Blanton M. R., Buries S. M., 
Cool R. J., Eisenstein D. J., Smith M. S. M., Wong K. C, Zhu 
G., 2012, ApJ, 746, 90 

Aird J., Nandra K., Laird E. S., Georgakakis A., Ashby M. L. N., 
Barmby P., Coil A. L., Huang J.-S., Koekemoer A. M., Steidel 
C. C, WillmerC. N. A., 2010, MNRAS, 401, 2531 

Best P. N., Kauffmann G., Heckman T. M., Brinchmann J., Chariot 
S., Ivezic Z., White S. D. M., 2005, MNRAS, 362, 25 

Bolzonella M., et al., 2010, A&A, 524, A76+ 

Brusa M., et al., 2010, ApJ, 716, 348 

Bruzual G„ Chariot S., 2003, MNRAS, 344, 1000 

Calzetti D., Armus L., Bohlin R. C, Kinney A. L., Koornneef J., 
Storchi-Bergmann T, 2000, ApJ, 533, 682 

Cao X., 2009, MNRAS, 394, 207 

Capak P., et al., 2007, ApJS, 172, 99 

Cappelluti N., et al., 2009, A&A, 497, 635 

Chary R„ Elbaz D., 2001, ApJ, 556, 562 



© 0000 RAS, MNRAS 000, 000-000 



18 E. Lusso et al. 



Ciotti L., van Albada T. S., 2001, ApJ, 552, L13 

Conroy C, Gunn J. E., White M., 2009, ApJ, 699, 486 

Dale D. A., Helou G., 2002, ApJ, 576, 159 

Di Matteo T., Springel V., Hemquist L., 2005, Nature, 433, 604 

Dickey J. M., Lockman F. J., 1990, ARAA, 28, 215 

Elvis M., Wilkes B. J., McDowell J. C, Green R. R, Bechtold J., 

Willner S. P., Oey M. S., Polomski E., Cutri R., 1994, ApJS, 95, 

1 

Fabian A. C, 1999, MNRAS, 308, L39 
Ferrarese L„ Merritt D., 2000, ApJ, 539, L9 
Gavignaud I., et al., 2006, A&A, 457, 79 

Goulding A. D., Alexander D. M., Lehmer B. D., Mullaney J. R., 

2010, MNRAS, 406, 597 
Haardt F, Maraschi L., 1991, ApJ, 380, L51 
Haardt F, Maraschi L., 1993, ApJ, 413, 507 
Hao H., Elvis M., Civano F, Lawrence A., 2011, ApJ, 733, 108 
Hao H., et al., 2010, ApJ, 724, L59 
Haring N„ Rix H.-W., 2004, ApJ, 604, L89 
Hasinger G, et al., 2007, ApJS, 172, 29 

Hoaglin D. C, Mosteller F, Tukey J. W., 1983, Understanding 

robust and exploratory data anlysis 
Hopkins P. F, Richards G. T., Hernquist L., 2007, ApJ, 654, 731 
Ilbert O., et al., 2009, ApJ, 690, 1236 

Just D. W., Brandt W. N., Shemmer O., Steffen A. T., Schneider 
D. P., Chartas G., Garmire G. P., 2007, ApJ, 665, 1004 

Kalberla P. M. W., Burton W. B., Hartmann D., Arnal E. M., Ba- 
jaja E., Morras R., Poppel W. G. L., 2005, A&A, 440, 775 

Kauffmann G., et al., 2003, MNRAS, 346, 1055 

Kawaguchi T., Shimura T., Mineshige S., 2001, ApJ, 546, 966 

Kelly B. C, Bechtold J., Trump J. R., Vestergaard M., Siemigi- 
nowska A., 2008, ApJS, 176, 355 

Kollmeier J. A., Onken C. A., Kochanek C. S., Gould A., Wein- 
berg D. H., Dietrich M., Cool R., Dey A., Eisenstein D. J., Jan- 
nuzi B. T., Le Floc'h E., Stern D., 2006, ApJ, 648, 128 

Komatsu E., Dunkley J., Nolta M. R., Bennett C. L., Gold B., 
Hinshaw G., Jarosik N., Larson D., Limon M., Page L., Spergel 
D. N., Halpern M., Hill R. S., Kogut A., Meyer S. S., Tucker 
G. S., Weiland J. L., Wollack E., Wright E. L., 2009, ApJS, 180, 
330 

Kormendy J., Richstone D., 1995, ARAA, 33, 581 

La Franca F, et al., 2005, ApJ, 635, 864 

Lilly S. J., et al., 2007, ApJS, 172, 70 

Lusso E., et al., 2010, Ac&A, 512, A34 

Lusso E., et al., 201 1, A&A, 534, Al 10 

Lutz D., et al., 201 1, Ac&A, 532, A90 

Lynden-Bell D., 1969, Nature, 223, 690 

Magorrian J., Tremaine S., Richstone D., Bender R., Bower G., 
Dressier A., Faber S. M., Gebhardt K., Green R., Grillmair C, 
Kormendy J., Lauer T, 1998, AJ, 115, 2285 

Mainieri V., et al., 2007, ApJS, 172, 368 

Mainieri V., et al., 2010, Ac&A, 514, A85 

Marchese E., Delia Ceca R., Caccianiga A., Severgnini P., Corral 

A., Fanali R., 2012, Ac&A, 539, A48 
Marchesini D., van Dokkum P. G., Forster Schreiber N. M., Franx 

M., Labbe I., Wuyts S., 2009, ApJ, 701, 1765 
Marconi A., Hunt L. K., 2003, ApJ, 589, L21 
Marconi A., Risaliti G., Gilli R., Hunt L. K., Maiolino R., Salvati 

M., 2004, MNRAS, 351, 169 
McCracken H. J., et al., 2010, ApJ, 708, 202 
Merloni A., et al., 2010, ApJ, 708, 137 

Merloni A., Rudnick G., Di Matteo T, 2004, MNRAS, 354, L37 
Netzer H, Trakhtenbrot B., 2007, ApJ, 654, 754 



Peterson B. M., Ferrarese L., Gilbert K. M., Kaspi S., Malkan 
M. A., Maoz D., Merritt D., Netzer H, Onken C. A., Pogge 
R. W., Vestergaard M., Wandel A., 2004, ApJ, 613, 682 

Pozzi F, et al., 2007, Ac&A, 468, 603 

Pozzi R, et al., 2010, Ac&A, 517, A11 + 

Richards G. T, et al., 2003, AJ, 126, 1 131 

Richards G. T, et al., 2006, ApJS, 166 

Salpeter E. E., 1964, ApJ, 140, 796 

Salvato M., et al., 2009, ApJ, 690, 1250 

Salvato M., et al., 201 1, ArXiv e-prints 

Scarlata C, et al., 2007, ApJS, 172, 406 

Schlegel D. J., Finkbeiner D. P., Davis M., 1998, ApJ, 500, 525 
Schmidt M., 1968, ApJ, 151, 393 
Scoville N, et al., 2007, ApJS, 172, 150 
Shankar F, 2009, New AR, 53, 57 

Shankar F, Salucci P., Granato G. L., De Zotti G., Danese L., 

2004, MNRAS, 354, 1020 
Shankar F, Weinberg D. H, Miralda-Escude J., 2009, ApJ, 690, 

20 

Shankar F, Weinberg D. H, Miralda-Escude' J., 2011, ArXiv e- 
prints 

Silk J., Rees M. J., 1998, Ac&A, 331, LI 

Silva L., Maiolino R., Granato G. L., 2004, MNRAS, 355, 973 

Simmons B. D., et al., 2011, ApJ, 734, 121 

Simmons B. D., Urry C. M., 2008, ApJ, 683, 644 

Steffen A. T, Strateva I., Brandt W. N., Alexander D. M., Koeke- 

moer A. M., Lehmer B. D., Schneider D. P., Vignali C, 2006, 

AJ, 131, 2826 
Steinhardt C. L., Elvis M., 2010, MNRAS, 402, 2637 
Tananbaum H, Avni Y., Branduardi G., Elvis M., Fabbiano G., 

Feigelson E., Giacconi R., Henry J. P., Pye J. P., Soltan A., 

Zamorani G., 1979, ApJ, 234, L9 
Tremaine S., Gebhardt K., Bender R., Bower G., Dressier A., 

Faber S. M., Filippenko A. V., Green R., Grillmair C, Ho L. C, 

Kormendy J., Lauer T. R., Magorrian J., Pinkney J., Richstone 

D., 2002, ApJ, 574, 740 
Trump J. R., et al., 2009a, ApJ, 700, 49 
Trump J. R., et al., 2009b, ApJ, 696, 1195 
Trump J. R., et al., 2011, ApJ, 733, 60 
Vanden Berk D. E., et al., 2006, AJ, 131, 84 
Vasudevan R. V., Fabian A. C, 2007, MNRAS, 381, 1235 
Vasudevan R. V., Fabian A. C, 2009, MNRAS, 392, 1 124 
Vasudevan R. V, Fabian A. C, Gandhi P., Winter L. M., 

Mushotzky R. F, 2010, MNRAS, 402, 1081 
Vasudevan R. V, Mushotzky R. F, Winter L. M., Fabian A. C, 

2009, MNRAS, 399, 1553 
Vignali C, Brandt W. N, Schneider D. P., 2003, AJ, 125, 433 
Young M., Elvis M., Risaliti G., 2009, ApJS, 183, 17 
Zamorani G., Henry J. P., Maccacaro T, Tananbaum H, Soltan 

A., Avni Y, Liebert J., Stocke J., Strittmatter P. A., Weymann 

R. J., Smith M. G, Condon J. J., 1981, ApJ, 245, 357 
Zheng W., Kriss G. A., Telfer R. C, Grimes J. P., Davidsen A. F, 

1997, ApJ, 475, 469 



© 0000 RAS, MNRAS 000, 000-000 



Bolometric luminosities and Eddington ratios ofAGN in XMM-COSMOS 



Table 7. Average Log A w values in Mbh - z bins for Type-1 AGN. 

LogM£ H <LogM BH > b N= bs N? orr (Logi^d) <r <Logi Edd > 
[ergs -1 ] [ergs -1 ] mean median 



z < 1.2 


7.00-8.05 


7.75 


21 


96 


-1.31 


0.49 


-1.65 


8.05-8.48 


8.36 


21 


26 


-1.28 


0.38 


-1.32 


8.48-9.10 


8.64 


18 


19 


-1.38 


0.46 


-1.37 


1.2<z<2.3 


7.70-8.34 


8.08 


36 


79 


-0.73 


0.32 


-0.82 


8.34-8.62 


8.50 


36 


52 


-0.90 


0.33 


-0.96 


8.62-9.30 


8.80 


37 


45 


-1.06 


0.45 


-1.07 


[a] Black hole mass intervals. 

[b] Median black hole mass. 

[c] Number of observed objects in each bin. 

[d] Number of predicted objects in each bin employing the V/V max method. 


Table 8. Mean Log /l Edd values in Mbh - z bins for Type-2 AGN. 




Log M£ H 
[erg s" 1 ] 


(Log M BH > b 
[erg s _1 ] 


N c 


N d 


(Log /tridd) 
mean 


(T 


(Log i Ed d> 
median 


z< 1.2 


5.45-7.70 


7.39 


104 


1046 


-1.46 


0.60 


-1.53 


7.70-8.07 


7.89 


107 


384 


-1.85 


0.60 


-1.94 


8.07-9.300 


8.28 


106 


262 


-1.89 


0.51 


-1.96 


1.2<z<2.3 


6.60-8.11 


7.83 


45 


201 


-0.91 


0.56 


-1.06 


8.11-8.47 


8.28 


45 


132 


-1.34 


0.48 


-1.40 


8.47-10.3 


8.69 


45 


65 


-1.30 


0.39 


-1.35 



[a] Black hole mass intervals. 

[b] Median black hole mass. 

[c] Number of observed objects in each bin. 

[d] Number of predicted objects in each bin employing the V/V max method. 
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